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Abstract 



A general electrodynamic theory of a grating coupled two dimensional elec- 
tron system (2DES) is developed. The 2DES is treated quantum mechanically, 
the grating is considered as a periodic system of thin metal strips or as an 
array of quantum wires, and the interaction of collective (plasma) excitations 
in the system with electromagnetic field is treated within the classical elec- 
trodynamics. It is assumed that a dc current flows in the 2DES. We consider 
a propagation of an electromagnetic wave through the structure, and obtain 
analytic dependencies of the transmission, reflection, absorption and emission 
coefficients on the frequency of light, drift velocity of 2D electrons, and other 
physical and geometrical parameters of the system. If the drift velocity of 2D 
electrons exceeds a threshold value, a current-driven plasma instability is de- 
veloped in the system, and an incident far infrared radiation is amplified. We 
show that in the structure with a quantum wire grating the threshold velocity 
of the amplification can be essentially reduced, as compared to the commonly 
employed metal grating, down to experimentally achievable values. Physically 
this is due to a considerable enhancement of the grating coupler efficiency be- 
cause of the resonant interaction of plasma modes in the 2DES and in the 
grating. We show that tunable far infrared emitters, amplifiers and gener- 
ators can thus be created at realistic parameters of modern semiconductor 
heterostructures. 

PACS numbers: 78.66.-w, 73.20.Mf, 42.55.Px, 78.45. +h 
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I. INTRODUCTION 



The motion of a fast electron beam across a periodic metal structure results in the 
radiation of electromagnetic waves. This phenomenon, often referred to as the Smith-Purcell 
effect ,0 provides the basis for a number of vacuum devices, such as travelling wave tubes and 
backward wave tubes. In these devices electrons accelerated by an applied electric field up 
to a velocity fdr, move in vacuum across a periodic metal structure (a grating or a spirale), 
which leads to an amplification or generation of electromagnetic waves at the frequency 
/ ~ Vdr/a, where a is the grating period. The drift velocity tidr here is determined by the 
applied electric voltage, so that these devices are voltage tunable amplifiers and generators. 

The vacuum devices successfully operate in the radio and microwave range. A further 
enhancement of the operating frequency presents severe difficulties because of the mechanical 
instability of a freely standing in vacuum periodic structure of metal wires with a very small 
period. The operating frequency of vacuum devices cannot therefore be extended up to the 
far infrared (FIR) range. 

In the late 70s an active experimental research of plasma oscillations in two-dimensional 
(2D) electron systems (ES's) in Si MOSFETs (metal-oxide-semiconductor field-effect tran- 
sistors) and GaAs/Ala;Gai_xAs heterostructures has been started.!^ In a considerable part 
of experimental work (for a review see Ref. ^ a FIR transmission spectroscopy has been 
used for the detection of 2D plasmons (Fig.|l], where v^r = 0). In this technique, the 2D 
plasmons which are normally nonradiative modes,0 are coupled to electromagnetic radiation 
by a metal grating placed in the vicinity of the 2D layer. An incident electromagnetic wave 
with the intensity Iq and the electric field polarized perpendicular to the grating strips, 
passes through the structure in the direction perpendicular to the 2D layer [z direction), 
and the spectrum of the transmitted wave T{uj) is registered. Well defined resonances which 
correspond to an excitation of 2D plasmons with reciprocal lattice vectors = {27cm/a, 0) 
are observed in the transmission spectrum (here a is the grating period and m is integer). 
In experiments of this type the energy of the external electromagnetic wave is converted to 
the energy of the 2D plasmons. 

Figure |I| resembles the geometry of the tunable vacuum amplifiers and generators: the 
system consists in a conducting electron (2D) layer where electrons can move under the 
action of an applied electric field, and an adjacent grating. First attempts to observe the 
emission of light from the grating coupled 2DES have been made in 1980.111 In these FIR 
emission experiments (Fig.|l], where the drift velocity Vdr 7^ 0, but the intensity of the 
incident light Iq is zero) a strong dc current is passed through the 2D layer (in the x direction 
perpendicular to the grating strips), and the emitted electromagnetic radiation is registered. 
The grating period in these solid-state structures can be made smaller than 1 fim, and the 
typical frequency of 2D plasmons falls in the terahertz (FIR) range. A successful realization 
of the 2D plasmon emission experiment could lead to a creation of a tunable solid-state 
sourse of the FIR electromagnetic radiation. In spite of the strong appeal of this idea and a 
number of more recent experiment s,EHli the intensity of radiation from the 2DES remains to 
be very small, and successfully working solid-state devices based on the discussed principle 
are absent so far. 

The energy of the dc current passing through an electron system in the presence of the 
grating is converted to electromagnetic radiation by two steps. First, it is transformed to the 
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energy of plasma oscillations in the beam by means of a current-driven plasma instability^ 
Then the energy stored in the plasmon field is converted to electromagnetic radiation by 
means of the grating. The plasma instability is developing in the system when the drift 
velocity of electrons exceeds a threshold value Vth, estimated as the plasma frequency of 
electrons in the beam divided by a typical grating wavevector Gi = 27r/a [see Eq. (pSf) 



below, as well as Ref. |T^. In vacuum devices the plasma frequency in the electron beam 
is much smaller than that in the 2DES, and electrons can be accelerated up to velocities 
much higher than those achievable in solid-state structures. The threshold condition for 
amplification is thus more difficult to satisfy in solid-state structures than in their vacuum 
counterparts, therefore the first attempts to realize the emission of light from the grating 
coupled 2DES have failed. 

The aim of this paper is to develop a general theory of the transmission, amplification 
and emission of light in the structure "grating coupler - 2DES" , and to find realistic ways 
for a reduction of the threshold velocity of amplification. We consider a propagation of 
light through the structure "grating coupler - 2DES with a flowing currenf (Figure |I| with 
fdr 7^ and Jo 7^ 0), and calculate the transmission T(u;, t>dr), refiection R{uj, t>dr), absorption 
A{ijj,vat) and emission E{u},VAr) (at Jq = 0) coefficients as a function of the light frequency 

the drift velocity fdr, as well as other physical and geometrical parameters of the system. 
In the literature the problem of the transmission of light (wdr = 0, Jq 7^ 0) has been solved 
analytically in a perturbative approachlllil (the grating has been treated as an infinitely 
thin metal layer with a weakly modulated density), and numerically in a nonperturbative 
approach.0 An emission of light from the structure "metal grating - 2DES" (fdr 7^ 0, Jq = 0) 
has been considered by Kempa et al^ using a numerical nonperturbative approach. We solve 
the problem analytically using the nonperturbative technique recently proposed in Refs. 19;20|. 
One of the main results of our work is that the amplification of waves can be drastically 
increased, and the threshold velocity can be essentially reduced (down to experimentally 
achievable values) in structures "2DES - quantum wire grating" (contrary to commonly 
employed structures with metal gratings). The effect is due to the resonant interaction of 
2D plasmons with plasmons of the grating, which leads to a remarkable enhancement of the 
grating coupler efficiency, and finally to an improvement of device characteristics. 

The paper is organized as follows. In Section |I| we develop a general theory of the 
scattering of electromagnetic waves on the structure "grating coupler - thin conducting 
layer" . In Section |T| we apply the general formalism to an analysis of the FIR transmission, 
refiection and absorption spectra of the grating coupled 2DES without the dc electric current. 
In Section |^ we study an amplification of FIR radiation passing through the system with 
a fiowing current. In Section [V| we discuss an emission spectrum of the structure (t>dr 7^ 0, 
Jo = 0) and compare our approach with that of Ref. |18[ In Section we summarize our 
results and formulate particular recommendations for a designing tunable solid-state FIR 
amplifiers. 



II. SCATTERING OF LIGHT ON A GRATING COUPLED THIN CONDUCTING 

LAYER: GENERAL THEORY 

In this Section we develop a general theory of the scattering of light on a structure 
"grating coupler - thin conducting layer" (Fig.|^). The grating coupler is treated as an 
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infinitely thini conducting layer with an electron density 

Ni{x)6{z) = — ak)6{z) 



(1) 



placed in the plane z = 0. The continuous function ni(x) is assumed to be zero at |x| > W/2 
and an arbitrary nonzero function at |x| < W/2, where W is the width of the grating 
strips and a is the grating period. The conducting layer (2DES) is placed in the plane 
z = D and described by the frequency and wave-vector dependent conductivity cr2D(q, i^) 
[all quantities related to the grating (2DES) will be supplied by the index 1 or ID (2 or 2D)]. 
Electromagnetic wave is assumed to be incident upon the structure along the z axis with 
the electric vector being polarized in the x direction, perpendicular to the grating strips. 
The system is infinite in the y direction, and a background dielectric constant e is uniform 
in all the space. 

The total electric field E*°* satisfies the Maxwell equations. 



2 TP tot 



V X (V X E*°*) 



4:71 d r 



c2 dt 



}'^ix)6iz)+nx)S{z-D) 



2D I 



with scattering boundary conditions at z ±oo. We search for a solution in the form 



[r,z 



cxt 
G 



[z) + E^''(z) 



AG-r—iwt 



(2) 



(3) 



where G = Gm = (27rm/a, 0), and the incident (external) and the scattered (induced) 
electromagnetic waves are written as 



E-*(z) 



|e^-v^^/^ E:^^ = EoSg,o, 

V 



and E^'^iz) = (^E'^i{z),0, E'^^{z)) . The field E'^^{z) satisfies the equation 



dz"^ 



^"""h&iz) + j!p{z - D)], 



ue 



and has a solution 



E'^%{z)U<o = AoexpiKcz), 



(4) 



(5) 



(6a) 



E'xGi^)\o<z<D = Bg sinh(KG2:) + Cq cofih.{KGz), 



where 



E':Uz)l^u = DGeM-^Gz), 



(6b) 
(6c) 

(7) 
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If Gm = 271771/ a = 0, the value k,g=o = —i'jJ\fllc is imaginary (the radiative boundary 
conditions at 2; ^ ±00 imply that Im kc* < 0), and the values Ag=o and Eq + -Dg=o give 
the amplitudes of normally reflected and transmitted waves. If kq is imaginary for several 
nonzero m (i.e. at oj^/e/c > Gm), the values and describe the amplitudes of reflected 
and transmitted waves in corresponding (m-th) diffraction orders. For all \m\ > a/A, where 
A = ^-Kcjuj-Je is the wavelength of light, Aq, and give the amplitudes of evanescent 
(nonpropagating) electric fleld. 

Using boundary conditions at the planes z = and z = D we relate the amplitudes of 
the electric fleld Aq, . . . , Dg to the Fourier components of the electric current: 



Ar 



Cr 



-kgD) 



(8) 



Be 



jl^G-Coexpi-KGD) 



(9) 
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G 
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cue 



jIg+J^g^M^^gD] 



(10) 



Together with the relation 



(11) 



between the current in the 2DES and the total electric fleld at the plane z = D we have 
flve equations for six unknowns Aq, . . . , Dg, j^^, J^'g- Using these equations we relate the 
total fleld at the plane z = to the current at the same plane. 



2'KiKG 

3x,G 



where 



(12) 



1 



e2D{G,uj) 



-2kcD 



(13) 



and 



2'KiKG ir-, 



(14) 



is the (relative) "dielectric permittivity" of the 2DES. 

Properties of the grating should now be introduced into the theory. Usually&^l one 
assumes the local Ohm's law for the grating, y^ix) = aiD{x,u!)El°^{x, z = 0), where the 
conductivity aiD{x,u!) is proportional to the local electron density (p. Then Eq. (|1^) is 
rewritten in the form of an integral equation 



Cf)2 uj'^f\ r fix' 



(15) 
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where = El°^{x,z = 0), -^{x) = ni{x) / {ni{x)) is a normalized electron density in a 

grating strip, the kernel L{x — x') is defined as 

^ ^ ooe ^ Kg ^ ^ 

f = W/a is the geometrical "filling factor" of the grating, and the angular brackets mean 
the average over the area of a grating strip, (...) = J{. . .)dx/W. 

A general scheme of solving Eq. (^) is presented in Appendix ^ Here we solve this 
equation approximately,ll2h0 assuming that the total (and induced) electric field inside the 
strips is uniform, = £'*°*(|x| < W/2,z = 0) = const. This approximation works very 

well in a metal grating if the frequency of electromagnetic wave is small as compared to the 
plasma frequency of the metal, and the width of the grating strips is large as compared to the 
Thomas-Fermi screening length. Under these conditions the electric field inside the strips is 
completely screened and i^^^'^"^ = const = 0. This approximation is also valid in a quantum 
wire (quantum dot) grating at an arbitrary frequency, if the wires (dots) are considered in 
an oblate cylindeiil (oblate spheroicS) model. This follows from the well known fact that 
an internal electric field in an arbitrary ellipsoid is uniform if the external one is uniform.!! 
This is also valid for wires or dots formed by a parabolic confining potential.i§ The validity 
of the model has been recently checked using a number of numerical approaches in Ref. 
P7| . It has been shown that the model gives reliable results for experimentally measured 
(macroscopic) values like for instance the transmission coefficient. 

Assuming that El^^^'^'^ = const, we get a relation between the total electric field inside the 
grating strips and the external field (see Appendix It has a form of a response equation 

TTiinside /l'7\ 

"cm- 

where the response function C(^^) is given by 



The form-factor in Eq. (|18|) 

«(G) = |(^?(a;)e'^-OI' (19) 
is determined by Fourier components of the equilibrium electron density in the grating strips. 



The response equation ([TtI) and the response function (|1^) are the main points of our 
theory. Having derived these equations we can now calculate fields and currents in all the 
space. In particular, for Fourier amplitudes of electric fields at 2; < and z > D which 
describe the refiected and transmitted field in all diffraction orders as well as the evanescent 
field we get 

-TT = -^G,o + M/ (Cr, ^) dG,o -r-. , (20) 



Dg _ , ^ 1 (, 2mfKG{aiD{uj)){d{x)e-'^-^) 
— -"G,o H 77;; r f G,o -r^. 



(21) 
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Equations (pO]), (pT|) give the general solution of the formulated problem. They have been 
derived under quite general assumptions and include both the electrodynamics of the grating 
coupler and nonlocal and quantum-mechanical effects in the response of the 2DES which 
enter the theory via an appropriate model of the conductivity a2D{G,uj). In subsequent 
Sections we apply the general theory to the problem of FIR response of the system "grating 
coupler - 2DES with and without a flowing current" . 



III. THEORY OF THE GRATING COUPLER: TRANSMISSION OF FIR 

RADIATION 

A. Approximations and preliminary notes 

Before applying the general results of Section H to the problem of FIR response of the 
grating coupled 2DES we specify the conditions of a typical experiment and make necessary 
approximations. First, we assume that the grating period a is small as compared to the 
wavelength of light A (in a typical experimental situation A ~ 300 fim, a ^ l/im). Under 
this condition only the G = components of the electric field describe outgoing waves, 
while all components with G 7^ are evanescent. The reflection r{uj) and transmission t{uj) 
amplitudes are then determined by the coefficients v4g=o and -Dg=o respectively, and we 
have 

.(.)^-f...(0,.)(f-^«)), (22) 



where €20(0, cu) = 1 + 2TTa2D{0,uj)/cy/e. The reflection, transmission and absorption coeffi- 
cients are then determined, as usual, by the relationsil 

Riu) = \r{u;)\^ T{u;) = \t{uj)\^ (24) 



A{uj) = 1- R{uj)-T{uj). (25) 

Second, we assume that the distance D between the 2DES and the grating is also small 
as compared to A. Then W{0,uj) = 62^(0,0;), and 

r{uj) = -l + t{uj), (26) 

where t^uj) is given by Eq. p3|). 

Third, we specify the model for the conductivity of the 2DES. We postpone an analysis of 
the nonlocal and quantum-mechanical effects in the 2DES to a subsequent publication, and 
describe the properties of the 2DES in the hydrodynamic modellH of cr2_D(q, c^;). Linearizing 
the continuity and Euler's equations^ for the density n and the velocity v of 2D electrons. 



7 



_ + V(nv) = 0, 



(27) 



9v 



+ (v- V)v 



-E - 72(V - Vdr) 



dt " ' m2 
where E = Eq + (5E, n = n2 + Sn, v = v^r + Sv, we get vjr = — (e/m272)Eo, and 



(28) 



m2 - qvdi.)(a; - qvj, + 272) ' 
For the average conductivity of the grating we assume, similarly, 



(29) 



nie 



nil uj + ill 



(30) 



In Eqs. (|28D-(|30D n^, rrij, and 7^ are the average electron density, the electron effective 
mass, and the momentum relaxation rate in the grating {i = 1) and in the 2DES {i = 2), 
respectively. 

Substituting the model expressions (^), ( ^0|) for the conductivity of the 2DES and the 
grating into Eq. (p3|) we get the following result for the transmission amplitudeEil 



00 + ii2 



1 



to + i'y2 + 'i^2\ tu + i'ji ({cu) J ' 
where the response function assumes the form 



X 1 



Here 



to + ill 

^li 
uj{uj + Z7i) 



1 + 



to + 272 
27r/nie^ 



E 



\GmHGm)e~^\^'-\''to%{G^j 



mieu^i (a; - G^t;dr)(cu - G^Vdr + ^72) - ^p2(G'm; 



tOp2{Gr^ 



'2TTn2e'^\Gr, 
m2e 



1/2 



is the frequency of 2D plasmons, Gm = 27rm/a, 

27r/nie^ 



mic-v/e 



(31) 



(32) 



(33) 



(34) 



27rn2e^ 
m2CA/e' 



(35) 



and 



8 



~,2 



m=f=0 



J2 \GrnHG^). (36) 



The physical meaning of the values Fi, and oOpi will be discussed in the Section [lllB 



Finally, we specify the density profile function 'd{x) which determines the form-factor 
a(G). In the main part of the paper we will use a semielliptic density profile, 

^ix) = ^[l-i2x/Wyf', (37) 

for which the form-factor a{Gm) is given by 

a{GJ = [2J,{z)/z]\ z = GmW/2, (38) 

where Ji in the Bessel function. In some cases we will also consider a step-like profile, 
^step(a:) = e{W/2 - for which 

"stcplGJ = (sinz/^)2, z = GrnW/2. (39) 



Using equations (|3^), as well as (p4D and ( P5|) one can show that the functions 

T{uj), R{uj) and ^(ci;) have one or more resonant features related with an excitation of 
plasma modes in the system. The resonance frequencies and linewidths depend on the drift 
velocity of 2D electrons fdr, as well as on other physical and geometrical parameters of the 
structure. If a resonant feature is well separated from others, i.e. when its linewidth is small 
as compared to the distance between adjacent resonances, the functions T^u), R^u) and 
A{uj) assume the following general form 

where Q is the resonance frequency, 7 and F are the nonradiative and the radiative decay 
rates, respectively. The total linewidth of the resonance is thus determined by the sum of the 
radiative and the nonradiative decay rates, the resonant values of the transmission, reflection 
and absorption coefficients, Tres = T{Q), i?res = R{^), ^res = ^(^), are determined by the 
ratio 7/F, 

72 F^ 27F 

= (TTTT' ^'""^ = (7 + F)2' ^''^ = (7 + F)2- ^^^^ 



The resonant values (|41|) characterize the strength of the resonant features. Note that 
the reflection (absorption) resonant amplitude is negligible as compared to the absorption 
(reflection) amplitude if F <^ 7 (F 7). In the following, we specify the values of fl, 7 and 
F in different considered cases. 

The rest of this Section is devoted to an analysis of the system without flowing current. 



i.e. at fdr = 0. In Section IV we analyze the general formulas at a flnite drift velocity. 
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B. Two limiting cases 



We start our analysis from two simple limiting cases, of a 2DES without grating coupler 
and of a grating without the 2DES. 



1. 2DES without grating 

If the grating is absent, then Fi = 0, and the transmission, reflection and absorption 
coefficients assume the form ( ^D] ) with 

j] = o, 7 = 72, r = r2. (42) 

The nonradiative contribution to the linewidth 7 is determined by the momentum relaxation 
rate of 2D electrons 72 and is due to the Drude absorption in the 2DES. The physical meaning 
of the value Eq. (|35D, is the radiative decay of oscillating 2D electrons in the 2DES 
taken in isolation (without the grating). Indeed, if one electron is placed in an electric field 
i?o exp(— iwt), it oscillates with an amplitude 5x ~ eEo/muj'^. This creates an oscillating 
dipole moment d ~ e^Eo/mco'^, which produces a dipole radiationiil with the intensity I ~ 
uj^(P/c'. Dividing the radiated intensity / by the average energy of the oscillating dipole 
W ~ muj'^dx^ one gets the radiative decay of a single electron Fq ~ e^u'^/mc'. When a 
sheet of electrons with an area density Ug is placed in an oscillating electric field, and the 
inter-electron distance is small as compared to the wavelength of light A, all N ~ n^A^ 
electrons within the coherence area ~ A x A radiate in phase. The average energy should 
then be multiplied by a factor of A^, while the radiated intensity by a factor of A^^. The 
radiative decay of an electron sheet is then given by the product TqN ~ TqHsX^ ~ n^e^ jmc 
in agreement with the exact expression (p5|). 



2. Grating without 2DES 



If the 2DES is absent, then To 



UJ, 



P2 



{Gr, 



0, and the transmission, refiection and 



absorption coefficients assume the form (^0]) with 



^ = ^pi, 7 = 7i, r = El. 



(43) 



The nonradiative contribution to the linewidth 71 is now due to the Drude absorption in 
the grating strips. The value Fi is proportional to the average electron density in the 
grating /ni, Eq. (|34D, and is the radiative decay of plasma oscillations in the grating 
taken in isolation (with the removed 2DES). The value Copi gives the resonance frequency 
of plasmons in a periodic array of grating strips (or quantum wires). Eq. (^) provides 
a functional dependence of ujpi on the equilibrium electron density d{r) in wires (similar 
functional dependencies for arrays of quantum dots and antidots have been found in Refs. 
|T^ and respectively). If the equilibrium electron density in strips has a semielliptic form 



Eq. (|36D gives (see Appendix | 



UJ. 



il = ^pl/5(/) = ^; 



pi 



24 



960 



(44) 
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where 



is the resonance frequency of plasmons in a single wire, and the factor /?(/) is due to the 
inter-wire interaction. If the wire is formed by an external parabolic confining potential 
Vext{x) = Kx'^/2, Eq. (^) reproduces an exact result Up-^ = K/rrii of the generalized Kohn 
theorem,!! see Appendix ^ 

Figure ^ shows the frequency dependencies of the transmission, reflection and absorption 
coefficients of a quantum wire array at two ratios of the collisional damping to the radiative 
decay, 71/ri ^ 1 (Fig. 0a) and 71/ri <^ 1 (Fig. ^b). In the former case the reflection of 
waves is negligibly small, and the transmission minimum is due to a peak in the absorption 
coefficient. In the latter case, the absorption of waves is small as compared to their reflection, 
and the transmission minimum is mainly due to the reflection peak. The width of the 
resonance in the second case is smaller than that in the flrst case (but does not tend to zero) 
and is determined mainly by the radiative decay. 

C. Grating coupled 2DES 

Now we consider the transmission of FIR radiation through a coupled structure "grating 
- 2DES", under the condition when no current is flowing in the 2DES {vdr = 0). We consider 
two different cases: the case of a metal grating, when the plasma frequency in the grating Upi 
is much larger (several orders of magnitude in a typical experiment) than the 2D plasmon 
frequency u!p2{Gi), and the case of a quantum-wire grating, when the plasma frequencies in 
the grating and in the 2DES are of the same order of magnitude. 

1. Metal grating 

Figure |^ demonstrates the frequency dependent transmission coefficient of the structure 
"metal grating - 2DES" at different values of the geometrical filling factor of the grating /. 
Three characteristic features are seen in the Figure. First, the position of resonances which 
we denote as Qu^m) does not coincide with the frequencies of the 2D plasmons ujp2{Gm), 
Eq. (|33D , shown in Figure |] by triangles for four lowest harmonics m = 1, 2, 3, 4. The index 
'12' here reminds that we deal with the coupled ID (grating) — 2D electron system. Second, 
the position and the amplitude of resonances for different m essentially depend on the filling 
factor / = W/a. At some values of / the amplitude of higher harmonics can be comparable 
with or even larger than those of lower harmonics (as is the case for the modes m = 1 
and m = 2 at W/a = 0.2, Figure or for the modes m = 2 and m = 3 at W/a = 0.6, 
Figure ^). At certain values of / some harmonics are not excited at all (e.g., the mode 
m = 2 at W/a = 0.6, Figure ^). Third, the amplitudes of the transmission resonances 
become smaller when the resonance positions fii2(m) approach the 2D plasmon frequencies 
^^p2{Gm), see, e.g., the evolution of the ^^12(1) mode amplitude with decreasing W/a (note 
the difference in the vertical axis scales in Figures ^a-c). 
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In order to understand these features we take the hmits Vdr = (no drift) and ni 



oo (metal grating) in the general formulas - (32). The resulting expressions for the 
transmission, reflection and absorption coefficients near the resonance u = Qi2{m) assume 
the form (|40|) where 

= QUm) = ujI,{GJ{1 - A J, (46) 

r = Tu{m) = ^ujl2{Gm)Am = r2-^^|m|A^, (47) 
and 7 = 72- The parameter 

= ^^H^iGm) exp(-2|G'„|D) (48) 

here depends on the harmonic index m and on geometrical parameters of the structure. 
Note that the physical parameters of the grating - the electron density rii, the momentum 
relaxation rate 71, and the effective mass mi - do not enter the formulas (P6| ) - (^HD, in 
which the grating is presented only via the geometrical parameters a, W, and D. The 
resonant values of the transmission, reflection, and absorption coefficients, Tj.es{m), Rresim), 
and Aj.es{m), are determined by Eq. (^I]) where 7 = 72 and T = ri2(m). 

As seen from Eqs. (^Bp, (0), and (|1T|), the resonance frequency, the radiative contribu- 
tion to the linewidth, and the strength of the resonance essentially depend on the parameter 
Am, which exponentially decreases with the distance D between the 2DES and the grating 
and oscillates as a function of the grating filling factor / [via the oscillating /-dependence 
of the form-factor a{Gm), Eqs. (|38|), (p9|)]. If Am tends to zero, the resonance frequen- 
cies f2i2(m) tend to the 2D plasmon frequencies ujp2{Gm), but the radiative decay and the 
strength of the resonances vanish, ri2(m) = 0, Ti.es(m) = 1. The value of A^ vanishes 
when vr/m coincides with any of zeros of the Bessel function Ji, for a semielliptic density 
profile (|37|) , or with any of zeros of the sine function, for a step-like profile. In order to get 
the maximum strength of the m-th plasma resonance one should thus satisfy the condition 
Ji (7r/m)=maximum (a semielliptic profile), or 

W W 

— m = 0.57, — m = 1.7, (49) 

a a 

etc. Figures ^ and |^ demonstrate the /-dependencies of the resonance frequencies Quim) 
[normalized by ujp2{Gi)], radiative contribution to the linewidth ri2(m) (normalized by r2), 
and the resonant transmission coefficient Tj.cs{m) for three lowest modes m = 1, 2, and 3 at 
D/a = 0.08, for the semielliptic density profile (p7|). The behaviour of Qi2{m), ri2(m) and 
^res("^) at the step-like profile is qualitatively the same. 



2. Quantum-wire grating 

If the grating coupler is made out of a metal, only the geometrical (but not the physical) 
parameters of the grating determine the observable transmission (refiection, absorption) res- 
onances. If the grating is made out of a 2D electron layer with similar plasma parameters 
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(the quantum-wire grating) the observable resonances are determined by both the 2D plas- 
mons in the 2DES, and the plasma modes in the grating. This gives additional possibilities 
to control the transmission spectra, especially in the finite drift velocity regime (Section |IV|). 

Figure 1^ shows the transmission coefficient T{u) of the structure ^^quantum-wire grating - 
2DES" at three different values of the grating plasma frequency ujpi. Geometrical parameters 
of the structure are the same as in Fig. |]c, where the transmission coefficient of the metal 
grating coupled 2DES is shown. Three new features are seen in Figure ^ as compared to 
Figure First, due to the presence of the grating plasmon resonance ujpi an additional 
resonance peak appears in the plot. Second, due to the interaction of 2D plasmons and the 
grating plasmon the resonance peaks are slightly shifted relative to their positions in Fig. 
Third, and the most important feature is a dramatic enhancement of the amplitudes of the 
2D plasmon resonances, in situations when the frequency ujpi approaches the 2D plasmon 
frequencies. Figure |]a,b. This effect is due to a resonant interaction of the grating plasmon 
with the 2D ones, and is especially pronounced for higher 2D plasmon modes, for which the 
amplitude of resonances is increased by about an order of magnitude (note the difference in 
the vertical axis scales in Figures ^ and |^). This effect is of a particular importance in the 
2DES with a ffowing current, as it allows one to increase the amplification and to reduce 
the threshold velocity in the structure with the quantum wire grating (Section |1V lj| ). 

IV. AMPLIFICATION OF WAVES 

Now we analyze the transmission of electromagnetic waves through a grating coupled 
2DES with a fiowing current. We start from the case of a structure "metal grating - 2DES". 

A. Metal grating 

Figure ^ shows the absorption coefficient of the structure "metal grating - 2DES" at 
relatively small values of the dc current {U = va^/vp2 < 1.6) in the frequency interval 
corresponding to an excitation of the first (m = 1) 2D plasmon harmonic. Here is the 
Fermi velocity of 2D electrons in the 2DES. The physical and geometrical parameters of the 
structure are the same as in Figure]^ {W/a = 0.6). The triangle at the bottom of the plot 
shows the position of the m = 1 2D plasmon harmonic (|33|). As seen from the Figure |^, in 
contrast to the case of the vanishing current {U = 0, lower curve), at finite drift velocities of 
2D electrons there are two modes, fli2{m, ±), associated with each harmonic number m. We 
label these modes by two indexes (m, ±), where the frequency of the + (— ) mode increases 
(decreases) with t>dr at small Vdr, dQi2{'m, +)/dvdT > and dQi2{m, —)/dvdT < at v^t 0. 
As seen from the Figure, 

lim ni2(m, +) = UJp2{Grn), (50) 

lim f2i2(m, -) = i7i2(m), (51) 

■Udr— >0 

where Qi2{m) is defined in Eq. (^61). The strength of the {171,+) mode vanishes when the 
drift velocity tends to zero. 
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In order to get a quantitative description of the resonant features shown on Figure |^ we 
take the hmit rii ^ oo (metal grating) in the general formulas (0) - (^). Assuming that 
I Co" + 272 1 ^ r2 and ujp2{Gm) ^ 72, and taking into account only the terms with m = ±|m| in 
the sum in Eq. (|32D , we find that near the (m, ±) resonance the transmission, reflection and 
absorption coefficients assume the form (|40|), where the resonance frequency Q = Qi2{m, ±) 
is determined by the equation 

nl,im, ±)= co^iGJ (1 - ^) + (GmVdrY 




the radiative decay T = Ti2{m, ±) is given by 

^ , ri2(m) [ Ujl^iGm)^^- i2GmVdr)^ 

2 I Up2{GJ^uj^p2iGm)Al + {AGmVa^ni - Aj2) 

and the nonradiative decay 7 = 72- Figure | shows the drift velocity dependencies of the 
resonance frequency f2i2(m, ±) and the normalized radiative decay ri2(m, ±)/ri2(m) for 
the mode m = 1 at parameters of Figure |^. The frequency fii2(m, +) [Qi2{m, — )] increases 
(decreases) with the square of the drift velocity at GmVdr ^ t^p2(Gm,)Am/4, and linearly, 

ni2{m, ±) ^ \uJp2{Gm)\/l - A^/2±GmVdr\, (54) 

at GmVdr ^ ^^p2 (Cm) Am/4 (the Doppler shifted plasma resonances). In the region 
^p2{Gm)V^ ~ Am < GmVdr < ^p2{Gm) the frcqucucy f2i2(m, — ) vanishes, and at GmVdv > 
^p2{Gm) it increases again with the positive slope, dfli2{m, —) / dvdr > at GmVdr > 
c<jp2(Gm)S The radiative decay, as well as the strength of the f2i2(m, +) resonance, equal 
zero at Vdr = and increase monotonously when the drift velocity increases (Figure ^). 
The radiative decay of the Qi2{m, — ) mode decreases from a finite value at Vdr = 0, vanishes 
at GmVdr = ^p2{Gm), and changes its sign at GmVdr > ^p2{Gm)- As seen from Eqs. (pj), 
when r equals zero, the reflection and absorption coefficients at the resonance uj = VL disap- 
pear, and the transmission coefficient equals unity. When F becomes negative, the resonant 
transmission coefficient exceeds unity, which means an amplification of waves, while the ab- 
sorption coefficient becomes negative, which means that the energy is transfered not from the 
electromagnetic wave to the electron system, but conversely, from the current driven elec- 
tron system (eventually from the battery which supplies the current) to the electromagnetic 
wave. The plasma mode f2i2(m, — ) thus becomes unstable at GmVdr > ^p2{Gm)^ Figure |^ 
demonstrates the drift velocity dependencies of the transmission, reflection and absorption 
coefficients at the resonance u = Qi2{m, ±) for m = 1. An amplification of the transmitted 
electromagnetic waves is explicitly demonstrated in Figure [1^ where we draw the frequency 
dependence of the transmission coefficient at different drift velocities at a larger range (as 
compared with Figure ^ of fdr (0 < t/ < 8). 

In the above calculations we have assumed that the resonance at = Qi2{m,±) is 
well separated from other resonances (i.e. that the width of the resonance line is small as 
compared to the distance to neighbour resonances). There are two effects in which this 
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approximation is insufficient. The ffist one concerns an accurate evaluation of the threshold 
velocity of the amplification of waves. We define the threshold velocity Uth by the condition 
Tres > 1. As follows from the above discussion and Figure ^ the resonant transmission 
coefficient exceeds unity when GmV dr > ^p2{Gm)- This inequality gives however only the 
lower estimate for the fth- In order to evaluate the threshold velocity more accurately one 
should take into account that the amplification of waves should exceed the Drude absorption 
in the 2DES (subsection [111 B 1| ) which is essential at low frequencies. Including this fact we 



get the following expression for the threshold velocity 



v,^ = ^^^^(l + ^X), (55) 



where the (positive) factor X is determined by the cubic equation 



872r2(272 + Ts) 
ri2(m)^2^(G„)A^ 



The threshold velocity thus consists in two contributions. The ffist one, 

can be reduced by chosing the structures with a low 2D electron density n2 and a small 
period a, and exploiting an excitation of higher 2D plasmon harmonics (here is the 
effective Bohr radius). The second contribution due to the correction AmX/2 in Eq. (|55[) 
has a complicated dependence on the density n2, momentum relaxation rate 72, the mode 
index m and the geometrical parameters of the structure. Qualitatively these dependencies 
can be understood if to note that the factor X equals A^^^, ii A ^ 1, and A^/^, ii A <^ 1, 
while the parameter A, in its turn, is proportional to 

/ 272 \ aexp{8GmD) 
A oc 72 1 + — a (58) 

(we consider the semielliptic density proffie). Thus, the second contribution to the threshold 
velocity can be reduced if the parameter W/a satisfies the conditions (^UD, the distance 
between the 2DES and the grating is small as compared to the width of the grating strips, 
D ^ W, and the grating period a, as well as the momentum relaxation rate 72, are taken 
to be as small as possible. The / dependence of the normalized threshold velocity Vth/vF2 
for several lowest modes m = 1, . . . , 4, and for parameters of Figure |] {ns2 = 3 x 10^^ cm~^, 
72 = 0.7 X 10^^ s~^, D = 60 nm) is shown in Figure 11. The divergencies of are related 
to zeros of the factor A^. 

The second effect which is not described by our single-resonance approximation is an 
anticrossing of modes with different m, which can be seen in Figure at f/ ^ 0.8 and 
at U ^ 1.5, where the relatively weak mode (3,—) intersects the modes (1,+) and (1, — ) 



respectively, as well as in Figure 10 at [/ ^ 5.0, where the mode (3, — ) (which has a 
positive slope with respect to v^r at so large U) intersects the mode (1,+) for a second 
time. The intersection points of modes (mi,—) and (m2, ±) (mi > m2) are determined 
by the relation f2i2(mi,— ) = f2i2(m2,±) [the modes f2i2(mi,+) and f2i2(m,2,±) do not 
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intersect at mi > 7712], and the transmission, reflection and absorption coefficients near the 
anticrossing can be found from Eqs. (|3T|) - (|3^) in the hmit ni — 00 (metal grating) if to 
take into account only the terms with m = mi and m = 1712 in the sum in Eq. (|32D. The 
most interesting situation is realized at a large drift velocity, v^r > ^^p2{Gmi)/Gmi, when 
an unstable plasma mode Qi2{mi, — ) intersects one of the stable plasma modes fli2{m2, ±) 
(mi > 1712). This occurs at the drift velocity 

(mi -),{m2,±) ^ ^p2(<Jmi) + ^^2(^*^2) 

and is accompanied by an enhancement of the amplification of waves, due to a resonant 
interaction of different plasma modes. In Figure |1^ one sees this effect at f/ ^ 5.0 [the 
intersection of modes (3, — ) and (1, +)]. A small resonance feature can be also seen at the 
intersection of modes (3, — ) and (1, — ) at f/ ~ 2.6 in the low-frequency range. 

In vacuum devices one can easily achieve the drift velocity sufficient for the amplification 
of electromagnetic waves. In solid-state structures ''metal grating - 2DES" the discussed 
values of the threshold velocity are rather large. In order to make realistic estimations of 
achievable drift velocities in semiconductor heterostructures with the 2D electron gas we 
refer to the paper of Wirner et a/.,!! in which the dependence of the average drift velocity of 
2D electrons as a function of the applied electric field has been experimentally investigated. 
The authors studied a GaAs/Al^-Gai-ajAs heterostructure with the density of 2D electrons of 
722 = 6 X lO^*' cm~^ and the (low-field) mobility of yU2 = 8 x 10^ cm^/Vs. The corresponding 
Fermi velocity of 2D electrons is fF2 = 1-06 x 10^ cm/s. The measured drift velocity 
of 2D electrons increases linearly with the applied electric field up to Eq ^ 50 V/cm, 
sublinearly at larger fields and then saturates at fdr ~ 1.8 x 10''' cm/s when the field is 
increased up to Eq 150 V/cm. Based on the results of this experiment we will assume 
that the really achievable experimental values of the ratio U = Vdr/vF2 in GaAs/Ala;Gai_2,.As 
heterostructures with the low-density high-mobility 2D electron gas are restricted by the 
value of U ^ 1.8. 

As seen from the above examples the threshold velocity is still well above the desired 
limit U ^ 1.8. It can be reduced, as compared to the numerical examples discussed above, 
by using smaller 2D electron gas density n2, smaller grating period a, and larger 2D plasmon 
harmonics m. One of the problems in using the higher m in the structures ''metal grating - 
2DES" is a small amplitude of 2D plasmon resonances with m > 1 and their rapid decrease 
with increasing m, see Figure ^. As we saw however in Section fll C 2| , the use of the quantum 
wire grating allows one to increase the amplitudes of the higher 2D plasmon resonances by 
an order of magnitude (compare Figures ^ and |^c), at the cost of the resonant interaction 
of 2D plasmons in the 2DES and the plasmons in wires. Using this effect, along with other 
methods discussed above, one can reduce the threshold velocity down to experimentally 
achiavable values. An amplification of electromagnetic waves in the structure "quantum 
wire grating - 2DES" is considered in the next Section. 



(59) 



B. Quantum-wire grating 

In order to make a realistic estimation of the transmission coefficient of electromagnetic 
waves in the structure "quantum wire grating - 2DES" we do this for a hypothetic sample 
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with parameters taken from published experimental papers. We assume that our sample is 
a GaAs/Al^Gai_^As heterostructure (mi = m2 = 0.067) with the density of 2D electrons 
in the 2DES of ^2 = 6 x 10^° cm~^ (taken from Ref. ^). The low- field mobility in Ref. p5 
was about /i2 ~ 8 x 10^ cm^/Vs which corresponds to 72 ~ 3.25 x lO^'' s~^. In the high-field 
regime {Eq ^ 150 V/cm) the mobility was by a factor of ~ 4 smaller, due to a heating 
of 2D electrons by a strong dc current. The dependence of the mobility on the dc current 
(or on the drift velocity) could be included into the theory through a phenomenological 
dependence of the momentum relaxation rate 72(7"e) on the electron temperature. For our 
estimations we use, for simplicity, the drift velocity independent Yalne 72 = 1.3 x 10^^ s~^, 
which corresponds to the mobility /i2 ~ 2 x 10^ cm^/Vs (roughly, this equals to the ratio 
Vdr/Eo at Eq ^ 150 V/cm in Ref. |35|) . Thus we assume the worst value of the momentum 
relaxation rate and take effectively into account the heating of 2D electrons by the strong 
dc current. For the grating we assume the same momentum relaxation rate of electrons, 
7i = 72- 

Choosing the geometrical parameters of the structure we have assumed that modern 
experimental technique allows one to create periodic microstructures with lateral dimensions 
of order of 0.1 fim, see e.g. Refs. |36|j37| . The width of the grating strips is therefore taken to 
he W = 0.1 /xm, while the period a = 0.175 fim is chosen in accordance with the rule ( ^Uj ) 
for m = 3. The distance between the 2DES and the quantum-wire grating is assumed to be 
D = 20 nm. 



Figures |T2| and demonstrate the calculated transmission coefficient of the structure 
"quantum wire grating - 2DES" , near the intersection point of the unstable 2DES plasma 
mode f2i2(3, — ) and the grating plasmon. Two different values of the electron density in the 



grating, ni = n2 = 6 x 10^° cm ^ (Figure 12) and rii = 2n2 = 1.2 x 10^^ cm ^ (Figure 13) 



are used. Three important features seen on Figures O and O should be mentioned. First, 
the resonant amplification of electromagnetic waves occurs at the drift velocities well below 
the experimentally achievable limit U ^ 1.8. The value of f/ ~ 1.4 (Figure |T2p corresponds 
to the drift velocity Vdr ~ 1.4 x 10^ cm/s and the dc current density jo ~ 0.13 A/cm in our 
example; in Ref. ^this velocity has been achieved at Eq ^ 50 V/cm. Second, the operating 
frequency of the amplifier lies in the vicinity of the intersection point of the grating plasmon 
(Dpi and the unstable mode Qi2{m, — ), at 



Upi + UJp2{Gm)\/'^ — ^m/'^ _ . . 

Vd, ^ , UJ ^ UJpi. (60) 

It is varied by the dc current (the drift velocity) within about 10 % with respect to ~ ujpi 
if the physical and geometrical parameters of the structure are kept constant (for instance. 



from ^ 0.73 to ~ 0.8 THz when U changes from 1.36 to 1.4 at Figure 12, or from ^ 1.03 to 



^ 1.13 THz when U changes from 1.52 to 1.58 at Figure |13D. The operating frequency can 
be also varied by changing the frequency Upi (compare Figures ^ and |T3|), in quantum- wire 
structures tunable, e.g., by a gate voltage. Third, the absolute value of the amplification 
of waves near the resonances can be as great as several tens of percents, which is due 
exclusively to the resonant interaction of the 2D plasmons with the grating plasmon in the 
quantum wire grating. For a comparison, in Figure |T^ we show the transmission coefficient 
of the structure ^^metal grating - 2DES" for the same parameters of the 2DES, the same 
geometrical parameters, and in the same frequency and drift velocity intervals as in Figure 
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T3| . A weak resonant feature which intersects the plot along the diagonal is the unstable 
mode f2i2(3, — ). As seen from Figure |14|, the amplification of waves in the metal grating 
structure is several orders of magnitude smaller than in the structure with the quantum wire 
grating (note the very large difference in the vertical axis scales in Figures |l^ and |l^. 

Thus an amplification of FIR radiation in structures with the quantum wire grating can 
be obtained at realistic, experimentally achievable parameters. 



V. EMISSION OF WAVES 

In previous Sections we have discussed an amplification of electromagnetic waves in the 
structure influenced by both the incident electromagnetic wave and the strong dc current 
(stimulated radiation). If the incident wave is absent, but the sample experiences a current 
flow, the system emits electromagnetic waves due to a disturbance of the thermal equilibrium 
(spontaneous radiation). This situation has been realized in experimental papers published 
so far. In order to describe the emission spectrum using the formalism developed above 
one should include into the consideration the equilibrium black-body radiation around the 
system, and take into account that the sample has a higher temperature (Te) than the 
environment (Tq). Taking into account that the sample reflects and transmits the incident 
black-body radiation with the intensity hbi^^, Tq), and emits the radiation with the intensity 
A{u!,Vdr)hb{^,Te)^ wc obtaiu that the emitted radiation registered by an external device 
(fllter) in the frequency interval between cu and cu + du is given byS 

E{uJ,Vdr) = A{uJ,Vd,)[IhbiuJ,Te) - Jbb(u;, Tq)] . (61) 

Here A{u!,Vdr) is the absorption coefficient of the structure calculated in Section Eq. 
(E5|), and 



is the intensity of the black-body radiation in the interval (a;, a; -|- du). 



Figure |T5| demonstrates the absorption (thin curves) and emission (thick curves) spectra 
of the structure "metal grating - 2DES" under the conditions of the experiment of Ref. ^ 
(n2 = 5.4 X 10^^ cm~^, a = 2 fim and a = 3 fim for two different samples, W/a = 0.6 in 
both cases, D = 62 nm, and the fllter linewidth A/ = 85 GHz). Plotting the Figure we 
have assumed a step-like electron density proflle in the grating, as the more relevant one in 
the case of the metal grating with wide strips, as well as the scattering rate 72 = 5 x 10^^ 
s~^, the environment temperature Tq = 4.2 K, and Tg = 100 K (Tg = 50 K) for the sample 
with a = 2 fim (a = 3 fim). The drift velocity was assumed to be vanishing. The whole 
behaviour of the emission spectra qualitatively agrees with the measured ones, the position 



of peaks in Figure |I5| is in a good quantatitive agreement with those measured in Ref. |T^, 
see Table |. 

It should be noted that the problem of the emission of light from the grating coupled 



2DES with a flowing current has been considered in Ref. |T8[ It has been solved in a complete 
analogy with the transmission problem. Such formulation of the emission problem is however 
not well deflned and cannot give the emission spectrum (|6lD (the fact that the emission of 
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light from the system is due to the black-body radiation of a sample heated by the dc current 
has been ignored in Ref. ^). Indeed, when the transmission of light is calculated, one gets 
a set of equations for Fourier components of the total electric field (Section |I| or Ref. 

E^CG'^G' =^G*, (63) 
G' 

where Mq,g' is an infinite matrix over reciprocal lattice vectors. The spectrum of eigen 
modes is determined by the equation det M = 0, the total self-consistent electric field can 
be found from Eq. (|63D , if the matrix M is inverted.0 The transmission problem is thus 
well defined. 

In Ref. 1^ the authors have derived a similar equation for the emission problem, when 
fdr 7^ and E°^* = 0. In this case however, the right-hand side of Eq. ( p^D vanishes, and 
the drift velocity v^r enters only the matrix M. The spectrum of eigen modes of the system 
(including unstable ones) can be calculated from the equation det M = 0, but the induced 
radiated field -Eg=o cannot be found in this fashion, as the right-hand side of Eq. ( |63D 
is zero. Instead of the intensity of the emitted waves, Kempa et al^ calculated the ratio 
of the macroscopic (radiated) field Egio to the microscopic (nonpropagating) field -Eq^o- 
This ratio characterizes the grating as a coupler of the plasmon field to the propagating 
electromagnetic radiation, but is not an appropriate characteristics of the emission process, 
as it does not vanish, for instance, at Vdi = 0. 



VI. SUMMARY 



We have developed a general analytic theory of the transmission, reflection, absorption 
and emission of electromagnetic waves in the structure "2DES (with and without the flow- 
ing dc current) - grating coupler" . We have demonstrated that an amplification, and hence 
a generation, of FIR radiation (the lasing effect) can be obtained in semiconductor mi- 
crostructures with realistic experimental parameters. Summarizing our results we formulate 
the requirements which should be met in order to create a successfully working amplifier of 
FIR radiation based on the grating coupled 2DES with a flowing current. 



1. The density of electrons in the 2DES should be small. This requirement seems to 
be paradoxical, as the intensity of the transmission resonance at the vanishing drift 
velocity v^r = becomes very small when n2 decreases. Nevertheless, at a large drift 
velocity the intensity of resonances becomes sufficiently large (Figure p!0|) , especially in 
the structures with a quantum wire grating (Figures |T^, 0), but the threshold velocity 



of amplification decreases with n2, Eqs. (|55D , (|57D . Note that this requirement has 



not been satisfied in previous emission experiments (for instance, in Ref. O the 2D 
electron density was by an order of magnitude larger than the value that we have used 
in Figures |T2|, |T^). Note also that in vacuum devices the plasma frequency in the 
electron beam used to be very small (Gifdr ^ i^p)- 

2. The mobility of 2D electrons should be sufficiently large, however this requirement 



is not so crucial as others. As we have seen in Figures ^ and |T3|, a considerable 
amplification of light could be achieved at a moderate mobility of /i2 ~ 2 x 10^ cm^/Vs. 
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The most important requirement imposed on the grating is that it must be a quantum 
wire grating, but not a commonly employediHli metal one. It is the resonant interac- 
tion of plasma modes in the 2DES and in the grating that allows one significantly to 
increase the amplification of light when using the higher 2D plasmon harmonics, and 
hence to reduce the threshold velocity of amplification. 

The grating period should be as small as possible. This condition seems to be in some 
contradiction with another one, a ^ D, which follows from the requirement that the 
interaction of plasma modes in the 2DES and in the grating [described by the exponent 
exp{—GmD), Eq- dH)] should be sufficiently large (for instance, in Ref. |l3|the grating 



period was by an order of magnitude larger than the one we have used in Figures |T^, 
[T3|) . Nevertheless, the inequality D <^ a (more accurately, D -C W , see below) should 
be considered as a condition for while the period should be taken to be small. As 
seen from Eqs. (|57D , (|58D , this leads to a reduction of the threshold velocity. 



The width of the quantum wires must satisfy the conditions (^) (or, Wm/a = 
1,2,3,..., if the steD-like profile seems to be more appropriate for a description of 
a particular systemEJ). As seen from Figure |TI|, a correct choice of the ratio W/a is of 
particular importance, especially when the higher 2D plasmon harmonics are used. 

The distance D between the 2DES and the grating should meet the condition D <C W. 
This requirement follows from maximizing the exponent exp(— Gm-D), Eq. (|48|), if to 
take into account that the ratio Wm/a is already fixed by the conditions (^91) . 



7. The operating frequency of amplifiers can be varied by the dc electric current fiowing 
in the 2DES and/or by changing the electron density rii (more generally, the resonance 
frequency Upi) in the quantum wire grating. 

We hope that these recommendations will finally help to create tunable FIR sourses and 
lasers, based on electron semiconductor microstructures of low dimensionality. 

ACKNOWLEDGMENTS 

I thank N. Savostianova for many useful discussions and W. L. Schaich for sending me 
the preprint of his paperS before publication. 



APPENDIX A: SOLVING THE INTEGRAL EQUATION (15) 



A general scheme of solving the integral equation (|T5D consists in the following. Let 
On{x) is a set of orthogonal polinomials with respect to the weight function 'dlx) which 
satisfy the condition 

dx 

—^{X)0n{x)0„^{x) = 5mn (Al) 

[for the profile (^) 0„(x) are the Chebyshev polinomials]. Substituting an expansion 



20 



E':\x) = Y.CnOn{x) (A2) 

n 

into Eq. (p!5[), multiplying by d{x)0„i{x) and integrating over dx we get an infinite set of 
equations 

+ Lmn)Cn = EoW{0, Uj)5mO, (A3) 

n 

which can be solved approximately truncating the matrix Smn + Lmn to a finite size N x N. 
Here 

L^n = ^^^^^^^^i^E^(G,^)«:G/3™(G)/3:(G), (A4) 
uje ^ 

(3m{G) = {^{x)Om{x)e^^'^) , and the star means the complex conjugate. The approximation 
accepted in the main body of the paper corresponds to = 1. Physically, truncating the 
matrix to the 1x1 size we neglect effects related to an excitation of quadrupole and higher 
eigen plasma modesH in the grating strips. 



APPENDIX B: RESONANCE FREQUENCY IN AN ARRAY OF WIRES 

Equation (|36D gives an expression for a dipole excitation resonance frequency in an array 
of wires as a functional of the equilibrium electron density i){x). It can be presented in 
different forms which clarify the role of the inter-wire interaction and the relation to the 
generalized Kohn theorem.!^ 



1. Inter- wire interaction 



In view of Eq. (^) one can write 



~,2 



pi 
'^1 
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\Gm\OL{Gm)- 



(Bl) 



For a semielliptic density profile (|37|) the form-factor a{G) is given by Eq. (pS]). Using the 
transformation 



E^(Gn.) = /$F(g)Ee--, 



(B2) 



where au = ak, we get 
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(B3) 



A;=— oo 



The term of the sum with A; = is independent of the grating period a and gives, after the 
integration over dq, unity. It describes the contribution of a single wire. The corrections 
due to the inter-wire interaction are then written as 
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l±) _±i rlT 

^ = 1 + 2^/ -J,2(/a:/2fc)(e- + e-"). (B4) 

OJpi Jo X 

In the first (second) integral in Eq. we rotate the integration path by an angle +7r/2 

(— 7r/2) to the upper (lower) complex half-plane. The function Ji is transformed to Ji, and 
we have 

The expansion (^4]) can now be easily obtained from Eq. (p5| ) at / ^ 1. 



2. Relation to the generalized Kohn theorem 

In many publications concerning the electromagnetic response of a confined system of 
electrons like quantum wires or quantum dots, the resonance frequencies are discussed in 
terms of not an equilibrium electron density, see Eq. (|36|), but of a confining potential. One 
can state a relation between these two approaches. 

Let an array of wires is formed by an external confining potential (potential energy) 

Ve^t{x) = 5I^ext(a; - ttk). (B6) 

k 

The total self-consistent potential Vtot(a;) = Vext{x) + Vind{x) is given by the sum of external 
and induced potentials, where the induced potential V^nd(a;) relates to the density (|T|) by the 
Poisson equation 

/ N Ane'^Ni(x) ^, , 
AVi,d(x, z) = -^S{z). (B7) 



Equation (BBf) can be written as 
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u;, = -— Y: \Gm\m,G^ f dxn,{x)e'''-' , (B8) 

^ mieWrii J cell 



where the integral is taken over an elementary cell. Due to Eq. (|B7| ) the Fourier component 
of the density A^i,g is related to the Fourier component of the induced potential V^nd,G = 
2ne'^Ni^G/^\G\, so that we can write Eq. (|B^ ) for a single wire in the form 

where A2 is the two-dimensional Laplacian and the integral is expanded onto the whole axis. 
Replacing Vind by the difference V^ot — Kxt, one sees that the contribution due to the total 
self-consistent potential vanishes, as V^tot is constant in points where electrons are. Thus, we 
have 

""^1 = ^^^;hv I ^^^i(^)^2v;xt(x). (BIO) 

In a parabolic confining potential Kxt(a;) = Kx'^/2, and Eq. ( [B10[ ) reproduces the statement 
of the generalized Kohn theorem,El 

col, = K/mi. (Bll) 
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FIGURES 



FIG. 1. The geometry of the considered structure. The system in infinite in y direction. 2D 
electrons are moving in the x direction perpendicular to the grating strips with the drift velocity f^r- 
A transmission spectroscopy experiment corresponds to Vdr = 0, Iq / 0, an emission spectroscopy 
experiment corresponds to Vdr ^ 0, Iq = 0, where Iq is the intensity of the incident wave. 

FIG. 2. The transmission, reflection and absorption coefficients of a quantum-wire grating 
with ni = 3 X 10^^ cm~^, rrii = 0.067 (GaAs), and two different values of the relaxation rate: (a) 
71 = 1.33 X IQii s-i (corresponds to 71/ri = 5.0), and (b) 71 = 0.53 x 10^° s^^ (ti/Fi = 0.2). 
The ratio W/a = 0.4. 

FIG. 3. The transmission coefficient of the structure "metal grating - 2DES" at three different 
values of the ratio W/a. Geometrical parameters: a = 1 /xm, D = 60 nm. Parameters of the 2D 
layer: n^2 = 3 x 10^^ cm'^, 72 = 0.7 x 10^^ s'^ m2 = 0.067 (mobility 1^2 = 375000 cmVVs). The 
grating parameters (n^i = 6 x 10^^ cm~^, 71 = 1.1 x 10^^ s~^, mi = 1) correspond to a typical 
(Au) grating coupler. The dielectric constant is e = 12.8. Triangles show the calculated positions 
of the 2D plasmon harmonics (^3|) for m = 1, . . . , 4. Note the difference of the vertical axis scales 
for different plots. 

FIG. 4. Normalized resonance frequencies i^i2('^)/'-^p2(G'i) for three different modes m = 1, 2, 
and 3 as a function of the geometrical filling factor / = W/a. The ratio D/a = 0.08. Note the 
difference of the vertical axis scales for different modes. 

FIG. 5. (a) Normalized radiative decay of the modes Qi2{m) and (b) the resonant transmission 
coefficient Ti.es(m), for three different modes m = 1, 2, and 3 as a function of the geometrical ffiling 
factor / = W/a. The ratio D/a = 0.08, the scattering rate of electrons in the 2DES (in Figure b) 
is 72 = 0.7 X 10^1 s-i. 

FIG. 6. The transmission coefficient of the structure "quantum wire grating - 2DES" at three 
different values of the grating plasmon frequency ujpi at a = 1 /xm, D = 60 nm, and W = 0.2 /xm. 
Parameters of the 2D layer and the dielectric constant e are the same as in Figure ^. Parameters 
of the grating: 71 = 0.7 x 10^^ s~^, m* = 0.067, the electron density: (a) Ugi = 4 x 10^^ cm~^, 
(b) Usi = 2.5 X 10^^ cm~^, (c) n^i = 1 x 10^^ cm~^. Triangles show the calculated positions of 
the 2D plasmon harmonics ( |33| ) for m = 1, . . . , 4, open triangles show the positions of the grating 
plasmon. 

FIG. 7. The absorption coefficient of the structure "metal grating - 2DES" at the frequency 
interval corresponding to the first (m = 1) 2D plasmon harmonic, and at small values of the 
dimensionless drift velocity U = VdT/vF2- Physical and geometrical parameters of the structure are 
the same as in Figure The black triangle at the bottom of the plot shows the position of the 
m = 1 2D plasmon harmonic (|33[). The weak mode which intersects the pronounced resonances at 
U w 0.8 and U w 1.5 is the mode (3, -). 
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FIG. 8. (a) The resonance frequency i}i2{m,±) and (b) the normahzed radiative decay 
ri2(?TT-, ib)/ri2(m) at m = 1 as a function of the dimensionless drift velocity of the 2D electrons 
U = Vdr/vF2- Physical and geometrical parameters of the structure are the same as in Figure 0. 



FIG. 9. The resonant values of (a) the transmission, (b) the reflection, and (c) the absorption 
coefficients for the modes Qi2{m,ziz) at m = 1 as a function of the dimensionless drift velocity of 
the 2D electrons U = Vdr/vF2- Physical and geometrical parameters of the structure are the same 
as in Figure 0. 

FIG. 10. The transmission coefficient of the structure "metal grating - 2DES" at the frequency 
interval corresponding to the first (m = 1) 2D plasmon harmonic in a wide range {0 < U < 8) of 
the dimensionless drift velocity U = Vdr/vF2- Physical and geometrical parameters of the structure 
are the same as in Figure ^a. 

FIG. 11. Threshold velocity of the amplification as a function of the grating filling factor 
/ = W/a for several lowest mode numbers m, and at = 3 X 10^^ cm-2, 72 = 0.7 x 10^^ s~S 
D = 60 nm. Thin lines show the first contribution to the threshold velocity (|57 



FIG. 12. Transmission coefficient of the structure "quantum wire grating - 2DES" for 
Til = n2 = 6 X 10^° cm~^, 71 = 72 = 1-3 x 10^^ s~^, a = 0.175 //m, W = 0.1 //m, and D = 20 nm. 
The frequency and the drift velocity intervals correspond to an intersection of the grating plasmon 
and the (3, — ) 2D plasma mode. The black triangle at the bottom of the plot shows the position 
of the grating plasmon (p6|). 



FIG. 13. The same as in Figure 12, but for ni = 2n2 = 1.2 x 10^^ cm ^. 



FIG. 14. The same as in Figure 13, but for a metal grating with Ugi = 6 x 10^^ cm ^, 
71 = 1.1 X 10^^ s~^, and mi = 1. Note the large difference in the vertical axis scale in this 



Figure and in Figure 13 



FIG. 15. The absorption (thin curves) and emission (thick curves) spectra of the structure 
"metal grating - 2DES" for parameters taken from Ref. |l^ The curves are vertically shifted for 
clarity, the absorption and emission are plotted in arbitrary units. 
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TABLES 



TABLE L Position of resonance peaks (THz) for two samples with the period a = 2 /zm and 



= 3 /xm, measured in Ref. 13 and calculated in this work, see Figure 15. For other parameters 
see the text. 
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a = 2 fim, theor. 
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a = 3 ^m, theor. 
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